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ABSTRACT 

We show that with the simple assumption of no correlation between the Lya equivalent 
width and the UV luminosity of a galaxy, the observed distribution of high redshift 
galaxies in an equivalent width - absolute UV magnitude plane can be reproduced. 
We further show that there is no dependence between Lya equivalent width and 
Lya luminosity in a sample of Lya emitters. The test was expanded to Lyman-break 
galaxies and again no dependence was found. Simultaneously, we show that a recently 
proposed lack of large equivalent width, UV bright galaxies (Ando et al. 2006) can be 
explained by a simple observational effect, based on too small survey volumes. 
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1 INTRODUCTION 

Two of the most common methods to find high redshift 
galaxies are based on detecting the Lyman-break or the 
Lyman-a emission-line of the galaxy. The Lyman-break 
technique, finding so-called Lyman-break galaxies (LBGs), 
has been extremely successful in gathering samples of galax- 
ies in the redshift range z ~ 3 — 6 (e.g. Steidel et al. 1996, 
1999, Pettini et al. 2001, Shapley et al. 2003, Bunker et al. 

2004, Ouchi et al. 2004). In parallel, the search for high red- 
shift galaxies by targeting the Lya emission of star-forming 
galaxies with narrow-band imaging, finding so-called Lya 
emitters (LAEs), has also been successful in the redshift 
range z « 2.5 - 7 (e.g. M0ller & Warren 1993, Cowie & Hu 
1998, Fynbo et al. 2002, Fynbo et al 2003, Matsuda et al. 

2005, Malhotra et al. 2005, Shimasaku et al. 2006, Tapken 
et al. 2006, Kashikawa et al. 2006, Venemans et al. 2007, 
Finkelstein et al. 2007, Nilsson et al. 2007, 2009, Ouchi et 
al. 2008). In both types of surveys, two of the most robustly 
measured properties are those of flux in the rest-frame ultra- 
violet (Muv) and the equivalent width of the Lya line (EW). 
For Lya emitters, these numbers are a direct product of the 
selection method. For Lyman-break galaxies, the Lya EW 
is only constrained in samples with spectroscopic follow-up. 

Lately, it has been suggested that there is a lack of high 
redshift galaxies with large Lya EWs and large UV fluxes. 
In Ando et al. (2006) a sample of Lyman-break galaxies 
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at z ~ 5 and z ~ 6 were studied for their Lya equiva- 
lent widths and rest-frame UV fluxes. Their results were 
also compared to Lya emitters at similar redshifts and they 
reported a lack of UV-bright objects with high EWs. The 
proposed explanations by the authors were those of different 
dust extinctions, smaller /larger amount of neutral hydrogen 
gas, age differences or gas kinematic effects. However, Ando 
et al. (2006) also suggested that larger samples are neces- 
sary to confirm this apparent deficiency of UV-bright, high 
Lya EW objects. Following this publication, several other 
authors have claimed to see the same effect. In the observa- 
tional publications of Ando et al. (2007) , Iwata et al. (2007) , 
Ouchi et al. (2008), Vanzella et al. (2009) and Shioya et al. 
(2009) a lack of these objects is reported. Ando et al. (2007) 
and Vanzella et al. (2009) both suggest that this may be 
a result of larger amounts of dust and a higher metallicity 
in UV-bright galaxies. Ouchi et al. (2008) propose that it 
infers a young, low-mass population of galaxies at high red- 
shift. The lack of high-EW, UV-bright galaxies has also been 
seen in theoretical work. Mao et al. (2007) present a simple 
physical model for high redshift galaxies, including a dust 
screen model. In parallel, Kobayashi et al. (2009) developed 
a semi-analytical model of high redshift galaxy evolution, 
including gas not homogeneously distributed, but rather in 
clumps. They both find their models to predict the lack of 
high-EW, UV-bright galaxies. On the other hand, there are 
also publications claiming not to see this effect. Verma et 
al. (2007) and Stanway et al. (2007) both study samples of 
z = 5 — 6 Lyman-break galaxies and see no UV dependence 



2 K.K. Nilsson, 0. Moller-Nilsson, P. M0ller, J.P.U. Fynbo & A.E. Shapley 



on the EW. Neither did Deharveng et al. (2008), observing 
nearly 100 Lya emitters at z ~ 0.3. Deharveng et al. (2008) 
also argue that the claim to see such an effect may be due 
to small number statistics. This is the same argument as Di- 
jkstra et al. (2007) propose. Their simple SFR-based galaxy 
evolution model does not reproduce any correlation between 
UV flux and Lya EW. 

In two recent publications, Pflamm-Altenburg et al. 
(2009) and Meurer et al. (2009) discuss how a varying IMF 
can cause a decline in Ha luminosity, and thus the EWs, at 
bright UV fluxes. Even though it is unclear what the rela- 
tion between Lya and Ha is in the high redshift Universe, 
there should exist a correlation between the two fluxes, and 
hence this effect could potentially also cause the Lya EW 
to decrease. Pflamm-Altenburg et al. (2009) predict that 
this effect will take place at star formation rates below 
0.01 M s yr" 1 . As surveys for high redshift galaxies are un- 
able to reach these low star formation rates, this particular 
effect is not expected to be observable. 

Given the large number of recent publications concerned 
with this issue, the number of proposed interpretations, and 
given that there still seems to be disagreement as to whether 
an effect is actually observed, we find it appropriate to sub- 
ject the issue to a rigorous statistical test. In this paper 
we address if there is a dependency between the Lya EW 
and UV luminosity (or Lya flux) in observed, high redshift 
galaxy samples. We focus on results from z ~ 3 Lya emitter, 
and Lyman-break galaxy, surveys as the amount of available, 
observational data is the largest at that redshift. Section [2] 
describes the method used. The results are given in sec. [3] 
and a discussion of the results is found in sec. [4] 

Throughout this paper, we assume a cosmology with 
Ho = 72 km s" 1 Mpc" 1 (Freedman et al. 2001), fi m = 0.3 
and f2 A = 0.7. 



2 METHOD 

The two observed parameters used here are the Lya and 
UV fluxes of a galaxy. For simplicity, we choose to perform 
our tests on plots involving the EW of the galaxy (see also 
Fig- [Dor Fig. [3} . This is a simple coordinate transformation, 
according to: 

EW Lya EE ^ (1) 
JUV 

where Fnya is the Lya flux and fuv is the flux density in 
the UV at the Lya wavelength. The 1-D EW distributions 
for high redshift galaxies, averaged over large samples, have 
been shown to closely resemble an exponentially declining 
function (Gronwall et al. 2007, Nilsson et al. 2009): 

P(EW)dEW = cons, x exp(-EW/w u )dEW (2) 

Here the w u is a scaling constant that can be fitted for for 
each data-set. The suggestion that there may be a lack of 
high luminosity objects with high EWs would then imply 
that Wu, at any given redshift, must be a function of the 
luminosity in the sense that the EW distribution becomes 
narrower for high luminosity objects as the high EW tail of 
the distribution is truncated. In principle this is easy to test 
simply by dividing a sample of galaxies into luminosity bins 



and then ask if there is any evidence for a significant change 
of w u between the high and low luminosity sub-samples. In 
practice the expected number of objects in the high lumi- 
nosity, high EW bin is too small to give a robust answer. 

Instead, we introduce a potential dependence on the 
EW from the Lya/UV flux, allowing a Monte Carlo ap- 
proach. To quantify the dependence, the parameter uew is 
introduced into the scaling constant according to: 

Wu.LAE = W0,LAE — UEW,LAE X j^jf^ (3) 

for Lya emitters, where wo,lae is the best fit scaling con- 
stant if uew,lae = and \^hya is the Lya luminosity in 
erg s _1 . For Lyman-break galaxies, instead we have: 

fuv 

W Ui LBG = W0, LEG ~ UEW,LBG X -^q-io 5 W 

where 

fuv = 10-°- 4 *< M ^+48.6) (5) 

The constants in Eq. [3] and [4] are carefully chosen to al- 
low a change in w u at the brightest end of the order ±Wo 
when —30 ^ uew ^ 30, and a change in the faintest end 
of the order a few A. These equations resolve into no de- 
pendence between Lya/UV flux and EW if uew = and 
a large dependence when uew — +~ ±Wo. A large positive 
uew will result in fewer objects with large EWs at bright 
fluxes, and a large negative uew in an excess of large EW 
objects at bright fluxes. Fig. [T] illustrates the effect of uew 
on the EW-flux distribution. It is important to note that 
it is necessary to fit the two types of galaxies with different 
equations, because they are flux limited in different parame- 
ters; Lya emitters in Lya flux and Lyman-break galaxies in 
the UV flux. Ideally, both galaxy samples would be fit with 
the same equation, with e.g. a dependency on only fuv- 
This would, however, infer a flux dependent lower limit to 
the EW of the Lya emitters, an effect which is clearly illus- 
trated in the left panel of Fig. [3] Because of this, we choose 
to fit Lya emitters and Lyman-break galaxies with separate 
equations and dependencies. 

To test the value of uew, we make a simulated dis- 
tribution of galaxies by drawing from 1-D distributions for 
the flux and the EW, with different values on uew, and 
test if the simulated galaxy distribution resembles the ob- 
served distribution. We do this in three steps. First we model 
the observed EW distribution, to find Wq. The procedure 
to find the best fit uew then includes creating simulated 
EW-flux distributions of galaxies by drawing randomly 
from the Lya/UV luminosity functions and the exponen- 
tially declining EW distribution function according to Eq. [2] 
and [3] or [4] for different values of Uew- The 1-D flux distribu- 
tions we draw from are the observed respective Lya and UV 
luminosity functions, summarised in Table [T] As this paper 
is focused on results at redshift three, the luminosity func- 
tion is drawn from Gronwall et al. (2007) in the case of Lya 
emitters. Several other luminosity functions have been pub- 
lished but they all agree very well (see also van Breukelen et 
al., 2005, Ouchi et al., 2008, Grove et al., 2009). In the case 
of Lyman-break galaxies we use the luminosity function of 
Reddy et al. (2008). 

Finally, we compare the simulated galaxies with the ob- 
served sample of 232 Lya emitters and 128 Lyman-break 
galaxies (see sec. [3] for a definition on the samples). For each 
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Figure 1. Plot of EW versus Lya luminosity for Lya emitters. Red stars mark the observed sample (see text for details). Each panel 
includes black contours with 68%, 95% and 99.7% significance of the simulated sample for different values of «ew,lae f° r LAEs. It is 
clearly seen that very large or very small u^w.LAE are ruled out by the observations. 



Table 1. Luminosity functions and EW distributions used in the creation of a random simulated sample. 



Type 


L* / M* 


<t>* 


a. 


Type of EW dist. 


■wo 




[ergs- 1 ] / [- 


[Mpc- 3 ] 






[A] 


Lya emitters 


42.66 


1.28 x 1CT 3 


-1.36 


Exp 


69 


Lyman-break galaxies 


-21.12 


1.12 x icr 3 


-1.85 


Exp 


29 



Luminosity function parameters are from Gronwall et al. (2007) for Lya emitters and from Reddy et al. 
(2008) for Lyman-break galaxies. Widths of equivalent width distributions (wo) are the one parameter fits 
to the total samples of Lya emitters and Lyman-break galaxies as described in the text. 
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simulated distribution, two statistical tests are performed 
to find the simulated distribution that best fits the observed 
data. Both methods are based on dividing the plane of data 
into four quadrants, with the intersection placed in an ar- 
bitrary point on the plane (see also Peacock 1983). In each 
quadrant the ratio between the number of galaxies in that 
quadrant to that of the number of galaxies in the total sam- 
ple is computed. This is done for both the data sub-set 
sample and the simulated galaxy sample. This process is 
re-iterated for every point in the plane until the largest dif- 
ference between the observed and test sample ratios is found 
in one quadrant. In our case the sampling of the plane has 
effectively one million resolution elements (1000 x 1000), en- 
suring that the maximum difference varies by less than 0.1% 
in the points nearest to the maximum. Based on the largest 
difference found, the following two tests are performed: 

(i) Monte-Carlo test 

(ii) 2-D Kolmogorov-Smirnov test 

In test (i) we merely find the largest difference between the 
simulated sample and the data. We calculate the errors on 
the ratios of test and observed samples and can thus deter- 
mine how dissimilar the two distributions are by calculating 
the difference in the ratios and the significance in the same. 
In test (ii) we perform the 2-D Kolmogorov-Smirnov test 
as described in Peacock (1983). The significance of the two 
distributions being similar is determined by: 



P(> Z 00 ) = 2x exp (-2(^ oc - 0.5) 2 ) 



/n x D n 



1 - 0.53n-°- 9 



niU2 
ni + m 



(6) 
(7) 
(8) 



where n\ and ri2 are the numbers of objects in the two re- 
spective samples and D n is the largest difference found in 
one quadrant. We end up with a set of simulated samples, 
one for each uew , each with a measure of how likely this 
sample resembles the real data. With these, we can find the 
distribution that best fits the data, and thus the best fit 
value of uew ■ 

Doing this once on the total sample of data will only 
reveal the best fit uew, without any understanding of the 
biases in the observed sample, or on the error bars. To un- 
derstand the full probability function of the uew parame- 
ter from the data we employ a jack-knife technique to the 
observed data. For each data-set (Lya emitters and Lyman- 
break galaxies) 3000 random sub-sets of data are created 
with the total number of galaxies in the total data-set mi- 
nus ten. For each of these data sub-sets, a new w u has to 
first be fitted for, after which the previously explained anal- 
ysis can be repeated. Thus, a best fit uew is determined for 
each of these 3000 data sub-sets, allowing an analysis of the 
distribution of uew- 



3 RESULTS 
3.1 Lya emitters 

To test the dependence of EW and Lya luminosity in Lya 
emitters, the observed samples of z ~ 3 Lya emitters of 
Gronwall et al. (2007), Nilsson et al. (2007), Venemans et 



al. (2007) and Grove et al. (2009) are used. A lower limit 
of restframe 20 A EW is set as this is typically the limit- 
ing EW for selection in narrow-band surveys. A lower limit 
to the Lya luminosity is also set to be equal to that of 
Gronwall et al. (2007); logE^ = 42.07 erg s _1 . After re- 
jecting objects in these catalogues with EWo < 20 A and 
log L^ya < 42.07 erg s 1 , and known AGN, we are left with 
a total sample of 232 LAEs with recorded Lya EW and lu- 
minosity. The uew,lae parameter is fit with values between 
—30 and 30 in steps of one. In Fig. [2]the histogram of best fit 
uew,lae to each of the 3000 random data sub-sets is shown. 
The median of the distribution is uew.lae ~ —4 and the 
mean is uew,lae = —3.66. Based on the skewed appearance 
of the histogram, the error bars are calculated by separately 
integrating the histogram to 68.3 % of each wing, split in 
the position of the median. The best fit uew,lae is then 



uew.lae 



1+9 



3.2 Lyman- break galaxies 

The largest sample of emission-line properties of Lyman- 
break galaxies has been published in Shapley et al. (2003). 
This is the sample that will be used here. An upper abso- 
lute UV magnitude is set to —20.0. EWs are constrained to 
be larger than 20 A in order to ensure spectroscopic com- 
pleteness. A total of 128 Lyman-break galaxies have been 
observed at redshifts 2.7 < z < 3.3 with EWo, Lya > 20 and 
Muv < —20.0. This population of galaxies is fit by select- 
ing an absolute UV magnitude from the luminosity func- 
tion of Reddy et al. (2008, see also Table [TJ, and then an 
EW from an exponential function, see Eq.[2]and Eq.[4] The 
uew,lbg parameter is fit with values between —30 and 30 
in steps of one. In Fig. [2]the histogram of best fit uew,lbg 
to each of the 3000 random data sub-sets is shown. The me- 
dian of the distribution is uew,lbg = 10 and the mean is 
uew,lbg = 10.23. The error bars are calculated in the same 
way as for Lya emitters, resulting in uew,lbg = lO^g . 



3.3 On the proposed lack of high-EW, UV bright 
high redshift galaxies 

The random Lya luminosities and EWs drawn for the sam- 
ple of test galaxies can also be converted to the restframe 
UV absolute magnitude, Muv, or vice versa. Thus, it is 
possible to reproduce the EW versus Muv plot (Ando et 
al. 2006) also for Lya emitters. In Fig. [3] we show this 
plot for the Lya emitter and Lyman-break galaxy samples 
and for uew = 0. The region of high-EW, UV bright ob- 
jects proposed by Ando et al. (2006) to lack objects was 
Muv < -21.5 (log /try > -10.84) and EWo > 20 A. As 
is seen in Fig. [3] this apparent lack of objects is easily ex- 
plained by small survey volumes. If there is indeed an inde- 
pendence between Lya luminosity and EW, then in order 
to find an object with high-EW and bright UV emission, a 
galaxy has to be drawn from the low probability ends of two 
distributions. The probability to find such an object is thus 
very small. In the example of uew = 0, and for galaxies with 
fluxes brighter than our simulation limits, the fraction of the 
total sample of galaxies in the quadrant with Muv < —21.5 
and EW > 20 A is 24 % for Lya emitters and 11 % for 
Lyman-break galaxies (with Muv < —20.0 and EW> 20). 
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Figure 3. Left Plot of EW versus fjjy for the Lya emitter sample. Black contours are 68%, 95% and 99.7% confidence contours of the 
simulated test sample, with uew,lae = 0, and the red points the observed points. The simulated sample is fully consistent with the 
observed sample. The lower limit to the distribution comes from the flux limitation in Lya. The upper left region of the plot is void of 
objects due to too small survey volumes. Right The same plot for Lyman-break galaxies, with contours from the simulated sample with 
v-EW,LBG = 0- 



The results of Ando et al. (2006) show two Lya emitters in 
this quadrant from a total of eight and three Lyman-break 
galaxies out of 17 in the total sample. This corresponds to 
percentages 25 ± 18% for Lya emitters and 18 ± 10% for 
Lyman-break galaxies, including very large error bars based 
on small number statistics, i.e. consistent with the simu- 
lated results. See also Table [2] for a summary on these re- 
sults. Based on the results for the sample from Shapley et 
al. (2003), there is a weak dependence between the EW and 
the UV flux, in the sense that there is a smaller fraction 
of high EW objects at brighter Mjjv , but uew,lbg = is 
only ruled out to Tier, or 73%. To confirm a real deficiency 
in objects in a quadrant would ideally require of the order 
100 objects in that quadrant, in this case requiring a to- 
tal sample of > 500 Lya emitters or > 1000 Lyman-break 
galaxies with confirmed fluxes and EWs. It is important to 
note that the analysis performed here was made at redshift 
z ~ 3 as most of the data exists at this redshift. Samples at 



Table 2. Percentages in high-EW, UV-bright corner in different 
samples. 



Type 


Simulations (uew = 0) 


Samples 


Ando 


LAE 


24% 


31±4% 


25±18% 


LBG 


11% 


9±3% 


18±10% 



Percentages of objects in the quadrant with Mjjv < —21.5 
and EWo > 20 A. Second column gives the percentages in the 
simulated samples with u^w = 0, third column in the observed 
samples used in this paper, and fourth column those of Ando 
et al. (2006). 

redshifts z ~ 5 — 6 are much smaller, and would produce a 
much weaker result. Several other factors make the higher 
redshift range more unsuitable for this test. Firstly, the mea- 
surement of the EW in very high redshift galaxies becomes 
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very uncertain due to the increasingly faint UV continuum, 
in some cases resulting in lower limits on the measured EW. 
Secondly, all line emission surveys, that are not spectroscop- 
ically confirmed, contain some fraction of low redshift inter- 
lopers. This fraction is expected to increase with increasing 
redshift (c.f. Kakazu et al. 2007). It is also uncertain if some 
neutral gas still remains at z ~ 6, which could affect the 
measurements of Lya EW/flux in unpredictable ways. For 
the sample studied here, at z ~ 3, these effects are not seri- 
ous. Until larger samples have been collected at either red- 
shift ranges, no conclusion can be drawn on whether there 
is a lack of high-EW, UV bright galaxies in the Universe. 



4 DISCUSSION 

How can the uew parameter be interpreted? A large posi- 
tive uew would indicate that there are fewer galaxies with 
large Lya EWs at brighter UV fluxes. A large negative uew 
indicates that it is more likely to find large-EW UV-bright 
galaxies. The next question is then what could cause this 
inflation or deflation of the EW distribution in the brighter 
UV flux slices? As Ando et al. (2006) argue, there could 
be several reasons for why there should be a dependence 
between these two measurable quantities. What is clear is 
that the time-scale on which the two quantities are sensitive 
to are different. The Lya flux, and thus EW, is only large 
in the first ~ 50 Myrs, assuming a single stellar population. 
The UV flux on the other hand can stay large for several 
hundred Myrs. In each galaxy, the UV and Lya flux is inte- 
grated over several star forming regions with different time- 
scales, with more massive galaxies having more star forming 
regions. If a larger fraction of these regions are in or out of 
their their Lya emitting phase, this could cause the EW dis- 
tribution to innate or deflate for UV-bright galaxies. Dust, 
metallicity, halo mass, infall rate, star formation rate, gas 
dynamics, and detailed gas morphology are also ingredients 
which could easily work to modify w u . 

In this paper we have searched for evidence of such 
an effect. We have done this in two ways. First we used 
a simple null-hypothesis of a single luminosity function and 
a constant w u for all luminosities. For a well defined sam- 
ple of 232 Lya emitters, the null-hypothesis predicts that 
56 of those should be in the (EW> 20 A, M uv < -21.5) 
section of parameter space. We find 72, which is 2a more 
than predicted by the simple model. This result is at vari- 
ance with earlier suggestions that this section of the EW- 
Muv plane is under-populated, if anything there might be 
marginal evidence that it is over-populated. Second, in or- 
der to parametrize this analysis in a more global way, we 
defined a description of the equivalent width distribution 
which allows it to change as a function of Lya luminosity 
via a parameter uew,lae- Following a Monte Carlo proce- 
dure we determined the best fit of a large number of random 
realizations to our sample and found uew,lae = -4+ 9 7 . This 
means that in a global sense the uew,lae parameter is con- 
sistent with zero for the entire sample, i.e. consistent with a 
constant w u . For the Lyman-break galaxies, the results are 
similar, with uew.lbg consistent with zero within 1.1c 

It is of interest to ask how strong the limits set by this 
result are. We will use a simple comparison to set those lim- 
its into context. In a recent paper (Nilsson et al. 2009) it was 



shown that w u is a function of redshift and that it has been 
found to drop from 69 A at z = 3 to 48.5 A at z = 2.25, i.e. 
a change of 20.5 A. From Eq.0we see that the most extreme 
effect, reducing the w u , occurs when uew.lae is as largely 
positive as possible. If we then include uew,lae = 5, the up- 
per envelope of the la probability for the Lya emitters, and 
compute what l^Lya corresponds to a shift in w u = 20.5 A, 
we find that this happens at a minimum Lya luminosity 
of 10 43 ' 11 erg s _1 , while 98.7% of the sample have smaller 
luminosities than this. In conclusion, our analysis has con- 
strained the dependence of w u on the Lya luminosity to be 
smaller than, or of the order of, the change observed as a 
function of redshift. 

The dependence- free case of uew = found here is 
intriguing, as it means that the Lya EW distribution is in- 
dependent on the Lya flux of Lya emitters and the UV flux 
of Lyman-break galaxies, and is then independent to all pa- 
rameters governing these fluxes (e.g. mass, age, metallicity, 
star formation rates etc.). In principle, the EW could still be 
correlated with the other respective parameter, i.e. the EW 
in Lya emitters could be dependent on Muv and vice versa, 
although the dependency between Lya flux and Muv m the 
galaxies would then have to be contrived. We consider such 
a scenario to be unlikely. 

The fact that our results are consistent with no depen- 
dence between Lya EW and the UV flux means that effects 
of varying the IMF, dust and age are not very important over 
the flux range that we probe in present day surveys of high 
redshift galaxies. That there is no age dependence indicates 
that the light observed from these galaxies is integrated over 
many star forming regions with varying ages. The analysis 
presented here would greatly improve by adding more sam- 
ples of both types of high redshift galaxies, but especially 
a larger sample of spectroscopically observed Lyman-break 
galaxies is necessary to finally draw any conclusions regard- 
ing the dependence between Lya EW and UV flux in high 
redshift galaxies. 
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